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Abstract — In this paper we consider the problem of transmit- 
ting a continuous alphabet discrete-time source over an AWGN 
channel. The design of good curves for this purpose relies on 
geometrical properties of spherical codes and projections of N- 
dimensional lattices. We propose a constructive scheme based on 
a set of curves on the surface of a 27V-dimensional sphere and 
present comparisons with some previous works. 

I. Introduction 

The problem of designing good codes for continuous al- 
phabet sources to be transmitted over a channel with power 
constraint can be viewed as the one of constructing curves 
in the Euclidean space of maximal length and such that its 
folds (or laps) are a good distance apart. When the channel 
noise is bellow a certain threshold, a bigger length essentially 
means a higher resolution when estimating the sent value. On 
the other hand, if the folds of the curve come too close, this 
threshold will be small and the mean squared error (mse) will 
be dominated by larger errors as consequence of decoding 
to the wrong fold. Explicit constructions of curves for analog 
source-channel coding were presented, for example, in [8| and 
0. 

In this work, we will consider spherical curves in the 2N- 
dimensional Euclidean space. We develop an extension of the 
construction presented in fl8) to a set of curves on layers 
of flat tori. Our approach explores geometrical properties of 
spherical codes and projections of iV-dimensional lattices. In 
the scheme presented here, homogeneity and low decoding 
complexity properties were preserved whereas the total length 
can be meaningfully increased. 

This paper is organized as follows. In Section II, we in- 
troduce some mathematical background used in our approach. 
In Section III we state the problem, while in Section IV we 
present a scheme to design piecewise homogeneous curves 
on the Euclidean sphere and describe the encoding/decoding 
process. In Section V we derive a scaled version of the 
Lifting Construction [6] suitable to our problem and present 
some examples and length comparisons with some previous 
constructions. 
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II. Background 

A. Flat Tori 

The unit sphere S 21 *' 1 c M. 2N can be foliated with flat tori 
H], Q as follows. For each unit vector c = (ci, c%, .., cjy) € 
R N ,d > 0, and u = (u u u 2 , . . . , u N ) E R N , let $ c : ->• 
R 2N be defined by 
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This periodic function $ c is a local isometry on its image, 
the torus T c , a flat iV-dimensional surface contained in the 
unit sphere S 2N ~ X . T c = & C (R N ) is also the image of the 
hyperbox: 

V c := {u E R N ; < u, < 27rc 4 }, 1 < i < N. (2) 

Note also that each vector of 5 f2W_1 belongs to one, and 
only one, of these flat tori if we consider also the degenerated 
case where some Ci is vanish. Thus we say that the family of 
flat tori T e and their degenerations, with c = (ci,C2, ..,cjv), 
||e|| = 1, ^ > 0, defined above is a foliation on the unit sphere 
of S 2 "- 1 C R 2N . 

It can be shown (Proposition 1 in [7|) that the minimum 
distance between two points, one in each flat torus TJ, and T c , 
is 

/ N 

d(T c ,T b ) = \\c-b\\ = [J2( 
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The distance between two points on the same torus T c given 
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by 



d($ c («),$ c («)) = 2 
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is bounded in terms of ||u — v\\ by the following proposition. 

Proposition 1. fi7§ Let c = (ci, C2, .., cn), ||c|| = 1, and let 
u,v € V c - Let A = ||u — v\\ and 8 = ||$ c (m) — 3> c (v)||. Then 
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(4) 
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where c^ = min Cj. 



B. Curves 

As pointed out in J51, some important properties of a curve 
from a communication point of view are its stretch and small- 
ball radius. Given a curve s : [a, b] i— > M. N , the Voronoi region 
V(x) of a point s(x) is the set of all points in R N which 
are closer to s(x) than to any other point of the curve. If 
H(x) denotes the hyperplane orthogonal to the curve at s(x), 
then the maximal small-ball radius of s is the largest r > 
such that B r (s(x)) n H(x) C V(x) for all x G [a,b], where 
B r (s(x)) is the Euclidean ball of radius r centered at s(x). 
Intuitively this means that the "tube" of radius r placed along 
the curve does not intersect itself. The stretch S(x) is the 
function ||s(x)|| where s(x) is the derivative of s(x). If ||s|| 
does not depend on x (as it is the case of the curves considered 
here) we will refer to the stretch as simply S. The length 
of a curve is given by J a S(x)dx. In this paper we will be 
interested in curves with large length and small-ball radius. 

III. Problem statement 

The underlying communication system we consider here is 
illustrated in Figure [T] Given an input real value x, within 
the unit interval [0, 1], the encoder maps x into a point s(x) 
of a curve s on S , which will be sent over an AWGN 
channel. By properly scaling the curve we guarantee that the 
transmitted energy is a 2 . The decoder will then compute an 
estimate for the sent value while trying to minimize the mean 
squared error (mse) E[(X — X) 2 ] of the process. 

For the torus layers scheme, the unit interval will be 
partitioned into M intervals of different length, and each of 
them mapped into a curve on one of the layers. It is worth 
noticing that, for the special case M = 1, if we choose the 
torus associated to the vector c = -4=(1, ■ . . ,1) to encode 
the information, then the scheme proposed here is exactly 
the one analysed in [8|. However, we are interested in the 
analysis for M > 1, in which case the curves to be presented 
outperform the proposed in 1 8 1 (see also [6|) in terms of the 
tradeoff between length and small-ball radius. 
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Fig. 1: Communication system 

The design of those curves in the next section is essentially 
divided in two parts. First, we choose a collection of tori on 
the surface of the Euclidean sphere S 21 ^^ 1 at least 6 apart. 
The approach for doing this is via discrete spherical codes 



in 



pN 



Second, we show a systematic way of constructing 
curves on each layer, via projection-lattices in l" -1 . At last, 
we give a description of the whole signal locus and summarize 
the encoding process. 



IV. Our approach 

A. Torus layers 

Given a fixed small ball radius 5, the first step of our 
approach is to define a collection T = {7i,T 2 , • • • ,?m} of 
flat tori on 5 2Ar_1 such that the minimum distance |3]l between 
any two of them is greater than 25. This step is equivalent 
to design a iV-dimensional spherical code SC with minimum 
distance 28 and consider just the points with non-negative 
coordinates. 

We denote this sub-code by 

SC+ = {c e SC, a > 0, 1 < i < N} . 

Each point c G SC+ defines a hyperbox V c ^ and hence a 
flat torus T c in the unit sphere 5 2Ar_1 . 

There are several ways of constructing spherical codes that 
can be employed here, e.g. 0, iU or even on layers of flat 
tori as introduced in Q. 

B. Curves on each torus 

Let T c be a torus represented by the vector c G SC+ as 
defined in the previous section. On the surface of T e we will 
consider curves of the form: 



(5) 



where C = diag(c 1; . . . , c£), u = uC = (ciUi, . . . , c^un), 
$ c is given by ([]} and x G [0, 1]. 

Provided that u G Z, N ,gcd(ui) = 1, those curves 
are closed (a type (m, . . . , Ujv)-knot), and due to 
periodicity and local isometry properties of $ c there 
lengths are 27r||«||. They are also the image through 
<£> c of the intersection between the set of lines 
W = {ux + n : n = nC,n G Z L } and the box V c - For 
c = -4=(1 , . . . , 1), these are exactly the curves analysed in 

m. 

Let r c (u) be the minimum distance between two different 
lines in W, then we have 
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where P^± (ft) denotes the orthogonal projection of ft onto 
the hyperplane u ± which is given by the standard projection 
formula 
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Let A c = C]Z © ... © catZ be the rectangular lattice 
generated by matrix C, then r c (u) is the length of shortest non- 
zero vector of the projection of A c onto u ± . Due to Equation 
Q, the small-ball radius 5 UtC of st c can be bounded in terms 
of r c (u) as follows: 
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where = mincj and c; > 0. Thus, for small values of 8 UC , 

we have <5 UjC « 7rr c (u). Our goal is to choose it in order to 
maximize r c (u). In addition, we also want to reach a contrary 
objective, which is the one of maximizing the arc length 
lu,c = 27T ||u|| of s Ta - 

It is possible to show (Proposition 1.2.9. in [5|) that the 
density of the lattice P^± (A c ), the projection of A c onto v,- 1 , 
is given by: 



A(P^(A C )) = 
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where An_i is the density of the best lattice in dimension 
(N — 1) and is the volume of the (N — 1) -dimensional 

unit sphere. For the case when all entries of c are equal, A c is 
equivalent to Z N and it was shown in [6 1 that we can make the 
above bound as tight as we want. We will show in Section V 
that this is also true for an arbitrary c i.e., that projections of 
the rectangular lattice A c can also yield dense lattice packings 
and therefore we can construct curves on the flat torus with 
the parameters arbitrary close to this bound. 



Example 1. Let N = 2. Consider the local isometry 
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on the flat torus T c and the line segment given by 
v(x) = xv, = x(2iruiCi, 2iru2C2),ui, u 2 G Z and < x < 1. 
The curve st c (x) will be the composition (&(v(x)) and we 
have 

C\C2 

y/w 2 cj + (W + l) 2 cl 

This curve in M 4 will turn around U\-times the circle 
obtained by its projection on the first two coordinates, whereas 
turning around u^-times the circle of radius c 2 given by its 
last two coordinates (a type {u\,U2) knot on the flat torus 
T c ). For this case, we can calculate the exact small-ball radius 
8 as d($(at)- L ), $(0)) where v 1 - = (—2ttc 2 U2,2ttciUi) and 
a = ircic 2 / \\v\\ 2 . In Figure [2] it is illustrated the curve 
(ui,u 2 ) = (4,5), with a > c 2 . 

C. Encoding 

Let T = {Tx,...,T M } be a collection of M tori as 
designed in Section |IV-A| For each one of these tori, let 
Sx k (x) = $(2iru k x), (k = 1,2, . . . , M) be the curve on T k , 
determined by the vector Uk §5§ and consider L = Ylj=x h< 
where l k is the length of Sr k - 

Now split the unit interval [0, 1] into M pieces according 
to the length of each curve: 



h = 
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for k = 1, 



Now, consider the bijective mapping 

fk : h [0, 1) 
fk(x) ~ - j 



x -E k j=l l i/ L 



Then the full encoding map s can be defined by 

s(x) := s Tk {fk(x)), if x e Ik- 



(12) 



The stretch of s will be constant and equal its total length 
L and the small-ball radius of s is the minimum small-ball 
radius S of the curves Sx k , provided that the distance between 
any pair of torus in T is at least 28. 

To encode a value x within [0, 1] we apply the map < fl2] >. 
The signal locus will be a set of M closed curves, each one 
lying on a torus layer Tj, and defined by a vector Uk- This 
whole process is illustrated in Figure [2] 
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Fig. 2: Encoding Process 

If the source is uniformly distributed over [0, 1], the en- 
coding process presented above is a proper one, since all 
subintervals will be equally stretched. For other applications, 
however, it could be worth considering another partition. 



D. Decoding 

Given a received vector y £ I 
decoding is finding x such that: 
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x = arg mm \\y 
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, the maximum likelihood 



s(x) 



Since exactly solving this problem is computationally ex- 
pensive we focus on a suboptimal decoder. 
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9, = arccos \^7^p» 1 <i < N. 

The process of finding the closest layer 
involves a iV-dimensional spherical decoding of 
7 = (71, 72, • • • , 7jv)> which has complexity 
0{MN). 

Let Ci = (cji,Cj2, • • • , Qjv) be the closest point in £ SC + 
to 7 and 

/ / 6i . &i \ ( n . N 

Vi = Cji cos — , sm — , . . . , c iN \ cos — , sm — 

V V 71 71/ V IN IN 

be the projection of y in the torus T Cj , i.e., 

< ||y-«||,V«eT Ci . 

From now on, we proceed the process by using a slight 
modification of the torus decoding algorithm [8| applied to 
the iV-dimensional hyperbox V Ci ■ The complexity of this 
algorithm is given by 0(N where Ui is the vector that 

determines the curve Hence, if M = ©(max^ Hujl^), the 
overall complexity of the process described in this section will 
be 0{N HujI^), the same as for the torus decoding. 

V. A Scaled Lifting Construction 

A. The construction 

The Lifting Construction was proposed in [6| as a solution 
to the problem of finding dense lattices which are equivalent 
to orthogonal projections of Z N along integer vectors ("fat- 
strut" problem). In this section we adapt that strategy in order 
to construct projections of the lattice A c which approximate 
any (N — 1) dimensional lattice (hence the densest one) with 
the objective of finding curves in approaching the bound 
Q. We adopt here the lattice terms as in ]2). 

The dual of a lattice A £ R N is the set A* = {x £ span(A) : 
(x, y) £ Z Vy £ A} where span(A) is the subspace spanned 
by a basis of A. Now let A c = ciZ © ... © cjyZ be the 
rectangular lattice generated by the diagonal matrix C (and 
with Gram matrix CC l = C 2 ). By scaling A c , we can assume 
that ci = 1. With this condition, if P u ± (A c )* denotes the dual 
of the projection of A c onto a vector u = (1, U2C2 ... , unCn) 
(v,i £ Z), then a generator matrix for P^±(A e )* is given by 
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what can be derived as a consequence of Prop. 1.3.4 J5) . 
Let now A be a target (N — 1) -dimensional lattice. Suppose 



that L* = is a lower triangular generator matrix for A*. 
We proceed to produce a sequence of lattices A^, w £ N, 
associated to the generator matrices 
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Similarly to |6|, it is possible to prove that L* w can be put 
on form ( p~3] > through elementary (integer) operations so that 
L* w is a generator matrix for P^± (A)* associated to a vector 
u depending on w. Also, we have (l/w 2 )L* w L*J; -> L* L* 1 as 
w —> 00 hence for large w, A will be approximated, in the 
sense of the Gram matrices, by projections of A c 

Example 2. Consider the hexagonal lattice R2§, which is the 
best packing in two dimensions and is equivalent to its dual 
and one of its generator matrix 



L* = 





V3 



We apply the construction above and reduce, through ele- 
mentary operations, the matrix L* w in order to put it on form 
( |13[ ). After re-scaling the rectangular lattice A c , we find the 
sequence of vectors 

u w = (c-l,-2wc 2 , {2w[wV3c 2 /ci\ - w)c 3 ). 

The projections of CiZ © C2Z © C3Z onto u ± will be, up to 
equivalence, arbitrarily close to A 2 . 

B. Comparisons: Curves in K 6 

Here we compare our approach of construct curves on torus 
layers with the curves constructed in [6] and [8] in terms of 
length for given small-ball radii. Given 8 > 0, we first consider 
a set of flat tori associated to a spherical code in R 3 , with 
minimum distance greater than 25, as described in Section 
IV-A Through the first inequality of ([8]), for each torus, we 
can find r c in order to guarantee that the curves on the flat 
tori will have small-ball radius at least 6 (this is also done in 
the case of the curves on the torus c = l/vZ(l, 1, 1)). We 
then look for the larger element of the sequence of vectors 
that produces a projection with minimum distance at least r c . 

In Figure [3] the first curve from the bottom to the top 
represents the exponential sequence J8). The second one is 
obtained by directly applying the Lifting Construction (6) to 
the hexagonal lattice and the last one displays the total length 
associated to our scheme. 

VI. Conclusion 

The problem of transmitting a continuous alphabet source 
over an AWGN channel was considered through an approach 
based on curves designed in layers of flat tori on the surface 
of a (2A^)-dimensional Euclidean sphere. This approach ex- 
plores connections with constructions of spherical codes and 




is related to the problem of finding dense projections of the 
lattice c{L © ... © c N Z. 

This work is a generalization of both the scheme proposed 
in ||8j and the Lifting Construction in J6). As a consequence, 
our scheme compares favorably to previous works in terms of 
the tradeoff between total length and small-ball radius, which 
is a proper figure of merit for this communication system. 

In spite of the improvements in terms of length versus small- 
ball radius, the constmctiveness, homogeneity and overall 
complexity of the decoding algorithm are features preserved 
from GO. 
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